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ABSTRACT 

Galaxy cluster merger statistics are an important component in understanding the formation of 
large-scale structure. Cluster mergers are also potential sources of systematic error in the mass 
calibration of upcoming cluster surveys. Unfortunately, it is difficult to study merger properties and 
evolution directly because the identification of cluster mergers in observations is problematic. We 
use large N-body simulations to study the statistical properties of massive halo mergers, specifically 
investigating the utility of close halo pairs as proxies for mergers. We examine the relationship between 
pairs and mergers for a wide range of merger timescales, halo masses, and redshifts (0 < z < 1). We 
also quantify the utility of pairs in measuring merger bias. While pairs at very small separations will 
reliably merge, these constitute a small fraction of the total merger population. Thus, pairs do not 
provide a reliable direct proxy to the total merger population. We do find an intriguing universality 
in the relation between close pairs and mergers, which in principle could allow for an estimate of the 
statistical merger rate from the pair fraction within a scaled separation, but including the effects of 
rcdshift space distortions strongly degrades this relation. We find similar behavior for galaxy-mass 
halos, making our results applicable to field galaxy mergers at high redshift. We investigate how 
the halo merger rate can be statistically described by the halo mass function via the merger kernel 
(coagulation), finding an interesting environmental dependence of merging: halos within the mass 
resolution of our simulations merge less efficiently in overdense environments. Specifically, halo pairs 
with separations less than a few /i _1 Mpc are more likely to merge in underdense environments; at 
larger separations, pairs are more likely to merge in overdense environments. 

Subject headings: cosmology:theory - methods:numerical - dark matter: merging histories - galaxies: 
clusters 



1. INTRODUCTION 

Galaxy clusters are of great interest in cosmology as 
they are the largest and most recently formed structures 
in the cosmological hierarchy. The clustering and num- 
ber density evolution of clusters are sensitive to both 
the growth function and the expansion history of the 
universe. Clusters contain a representative sample of 
baryons and dark matter and are thus also fascinating 
laboratories in which to study the influence of baryonic 
physics on the formation of large-scale structure, sub- 
structure, and its effect on the gravitational potential 
and dark matter halo shape. Statistical measures such 
as the mass function, the rate of structure growth, and 
the clustering of large structures, are fundamental pre- 
dictions of cosmological models. Cluster observations are 
therefore expected to provide some of the most important 
constrai nts on fundame ntal cosmology and astrophysics 
(see e.g. iBorganil (|2006) for a recent review). 

Cluster formation histories are frequently punc- 
tuate d by large jumps in mass from major mergers 
(e.g. iCohn fc White! J2005)). These mergers are one 
of the primary mechanisms for the buildup of mass in 
clusters and superclusters. In the standard paradigm 
(lKaiserlll98'llEfstathiou et al.lll988HCole fc Kaiserlll989t 
iMo fc White! 119961 : ISheth fc Tormenl 1 1999h . observable 
properties such as the degree of spatial clustering 
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depend only on the cluster mass. However, recent 
theoretical studies indicate that many cluster observ- 
ables, such as spatial clustering, concentration, galaxy 
velocity dispersion, gas distribution and its attendant 
observables such as X-ray emissions and SZ decrement, 
depend on the cluster's formation time, mass accretion 
history, large scale en vironment (collectively referred to 
as "as s embly bias") (IWechsler etall l2002t IZhao et al. 
2Q0M ISheth fc Tormenl 120041: lOao. Snringel fc Whitf 
2005L IWechsler et al.| [200l ICroton Gaq fc White [20061: 
Harker et alj 120061: IWetzel et al.l 120071 : iGao fc Whit 
IWang. Mo fc Jind 
Hahn et alj 
Sandvik et alJ 120071) . 



120071 Uing. Suto. fc Mol 
" 120071 : iMaccio et all [2003 



In addtion, there is a de 
pendence on recent merger history ( "merger bias" ) 
(IScannapieco fc Tha cker 2003; Ro wley. Thomas fc Kavl 
' 20041: iFurlanetto fc Kamionkowskil 120051: IWetzel et all 
20071: iPoole et aLl 120071: Ueltema et all 120071 : 
Hartley et all I2007D . In addition, recent studies 
have claimed observational dete ction of assembly - 
bias, though with mixed results (|Berlind et al.1 120061 : 
lYang. Mo. fc van den Boschl [200(1 . Theory can predict 
cluster properties as a function of their mass, which is 
dominated by dark matter and thus cannot be directly 
measured. Since all methods of observing clusters are 
sensitive to the effects of assembly/merger bias, it is 
necessary to develop a more detailed understanding 
of the mechanisms of structure formation. Specifying 
a correlation function and mass function, and the 
evolution of these quantities, may not be sufficient to 
connect theory to observation. Understanding cluster 
merger properties is therefore crucial to utilize these 
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objects as probes of cosmology. 

We focus prim arily on galaxy cluster mergers. 
IWetzel et al.l ([2007) used a large simulation volume to 
probe high-mass halos with good statistics. They found 
that halos of mass M > 5 x 10 13 /i _1 M that have un- 
dergone a recent (within 1 Gyr or less) major merger or 
large mass gain exhibit an enhancement in their spatial 
clustering of up to ~ 10-20% on scales of 5-25 /i _1 Mpc 
compared with the entire halo population at the same 
mass. They noted that this merger bias persists for the 
rcdshift range < z < 1, and that the bias increases with 
larger merger mass ratios or shorter merger timescales. 
If this merger bias remains unaccounted for, then mass- 
observable relations that connect theory with observation 
calibrated on the basis of clustering may be suspect. 

Moreover, the gas properties of a recently merged halo 
cluster can be quite different from the general popula- 
tion. This can be mitigated by selecting "relaxed" halos 
(presumably those that have not had recent major merg- 
ers) when calibrating observables. However, the scaling 
relations from this selected sample may be biased with 
respect to that from the overall cluster population. In 
addition, if a merger accidentally entered the "clean" 
sample (e.g. because it occurred along the line of sight), 
it could substantially bias the calibration, and hence the 
ensuing mass determination. It is thus important to pre- 
dict the fraction of halos which have had recent merger 
activity. Observing cluster mergers as a function of red- 
shift, studying the merger bias, and correlating mergers 
with other tracers of density and environment will also 
shed light on the nature of st ructure forma t ion. 

A nu mber of authors ( e.g. iPatton et all (119971 I2000L 
20021): IPePropris et alJ (l2005j)JInfante et all (|00| ): 
Lin et al.l (12004): jC onsclice 1120061) : feell et all (j200ff i: 
Masiedi et all (|2006D : iDe Propris et alJ ([20071 )) have 
studied the statistics of galaxy merging using close pairs 
as a proxy for mergers. In this work we extend those 
analyses to the cluster scale. It is important to quantify 
the co nditions for which the pair proxy assumption is 
valid. iBerrier et al.l (|2006[ ) used simulations to examine 
the utility of galaxy pairs in measuring the redshift evo- 
lution of halo merger rates out to z ~ 3. They tracked 
the formation and evolution "host" halos of dark mat- 
ter, as well as self-bound density peaks within host ha- 
los, i.e. "subhalos," a term which includes both satellite 
objects and the central host itself. Using the assump- 
tion that galaxies are found at the centers of subhalos, 
they found that galaxy pair counts can probe the rates 
at which satellite galaxies merge with the central galaxy, 
i.e. subhalo mergers. These pairs can also constrain the 
galaxy Halo Occupation Distribution, i.e. the statistical 
number of galaxies found within a host halo of a given 
mass. However, galaxy pairs cannot be used to measure 
the merger rates of separate massive host halos. 

The goal of this paper is to investigate whether 
close spatial pairs of halo objects (galaxy clusters) 
can be used to observationally study the properties of 
mergers. Identifying a cluster merger in an observa- 
tion is not straightforward because it relies principally 
on the object appearing mor phologically "disturbed" 
([Rowley. Thomas fc Kavll200l ). Furthermore, it is dif- 
ficult to estimate the completeness and the contamina- 
tion of such an observed merger population. If close pairs 
prove to be a sufficient proxy for mergers, they would pro- 



vide a clean tool to study merger statistics. We concen- 
trate primarily on massive halos (M > 5 x 10 13 /i _1 M Q ), 
corresponding to galaxy groups and clusters, in contrast 
to previous work which has focused on pairs as a proxy 
for galaxy mergers. However, our results also apply to 
any class of object that singly occupies a halo, particu- 
larly field galaxies at high redshift. 

In $2] we describe the simulations and halo catalogs 
used in our study. In [J3]we present statistics for merger 
rates as a function of halo pair masses, pair separations, 
and merger time interval to quantify the extent to which 
close pairs can be used as proxies for cluster mergers in 
the best-case scenario where halo spatial positions are 
known accurately. In S2]we consider observational com- 
plexities such as scatter in the halo mass, redshift space 
distortions, and redshift space errors. In [J5]we present a 
mass-dependent fit to the merger kernel and examine how 
large-scale density impacts the efficiency of pair merging. 
In f}6] we discuss the prospects for studying merger bias 
using close cluster pairs. We conclude in fj7l 

2. SIMULATIONS 

Our study is conducted using two high-resolution N- 
body simulations of cold dark matter in a flat ACDM 
cosmology with parameter values VIm = 0.3, FIb = 0.046, 
h — 0.7, n = 1 and os = 0.9. Our sim ulations employ the 
HOT code ([Warren fc Salmonl[l99l in a (1.1 h^Gpc) 3 
and a (2.2 /i _1 Gpc) 3 volume with periodic boundary con- 
ditions, using a Plummer softening length of 35 h~ x kpc. 
Gaussian initial conditions were randomly generated for 
the 1024 3 particles of mass 10 11 h~ 1 M Q (smaller simula- 
tion) and 8 x 10 11 h~ 1 M@ (larger simulation) at an initial 
redshift of z = 34. Simulation outputs were stored in in- 
tervals of 1 Gyr between redshifts z rj 1 and z = 0, 
with the last interval of each simulation being shorter: 
0.6 Gyr (smaller simulation) and 0.5 Gyr (larger simula- 
tion). All time intervals cited below represent the total 
time elapsed between two simulation outputs. Thus a 
merger timescale of At — 1 Gyr indicates that two sepa- 
rate halos have merged within 1 Gyr, an upper limit to 
the actual time to merger. 

We generate a halo catalog for eac h output using the 
Friends-of-Friends (FoF) algorithm ([Davis et alJ 1 19851) 
with a linking length of b = 0.15 of the mean inter- 
particle spacing. These groups correspond to a density 
threshold of ~ 3/(27r6 3 ) and enclose primarily virial- 
ized material. In this work we will quote FoF masses, 
which are about 20% smaller than "virial" masses (cor- 
responding roughly to FoF masses with b = 0.2, which 
are m ore commonly found in the literature; see White 
([200lh for more details). Since we are examining merg- 
ers, we use a smaller linking length to decrease contam- 
ination that arises from close neighboring halos being 
bridged by a narrow string of particles (although using 
a larger linking length of b = 0.2 changes our results by 
only a few percent, see Fig. 0J. The halo catalogs of 
the smaller simulation include all halos of mass greater 
than M k 5x 10 12 h~ l M Q (> 50 particles), though 
in our study we are concerned primarily with halos of 
mass M > 5 x 10 13 /i~ 1 Af Q , of which there are around 
75,000/i" 3 Gpc" 3 at z = 0. Our larger (2.2 h^Gpc) 3 
but less resolved simulation catalogs include all halos of 
mass greater than M « 4 x 10 13 h M@ (> 50 particles). 
We do not consider substructure within host halos. 



Close Pairs as Proxies for Mergers 



3 




Progenitor Separation (Mpc/h) 

Fig. 1. — Number density of progenitor pairs as a function of 
(comoving) binned separation. We consider children with masses 
above 5 X 10 13 h —1 Mq and restrict ourselves to cases where the 
two most mass-contributing progenitors have total mass ratios 
M2/M1 > 0.2 (major mergers). A fixed merger time interval of 
At = 1 Gyr is used for child halos at 2 = 0.04 (solid), z = 0.44 
(dotted), z = 0.75 (dot-dash). 

Merger trees were constructed from the set of halo cat- 
alogs by specifying a parent-child relationship, where a 
"parent" is any halo that contributed mass to a halo at a 
later time, i.e. a "child." We define a parent contributing 
more than half of its mass as a "progenitor." Under this 
restriction, a progenitor can never have more than one 
child. We define a merger as a child having more than 
one progenitor. In cases where we select on a child that 
had more than two progenitors, we apply the two body 
approximation by considering only the two progenitors 
that contributed the most mass to the child. Using in- 
stead the two most massive progenitors, as opposed to 
the two most mass-contributing progenitors, changes the 
progenitor identification in less than 1% of all mergers 
that we consider. In addition to this progenitor merger 
tree, a list of all contributing parents (not just progen- 
itors) is stored for each child. This flexible storage of 
parent data allows us to study two body mergers, merg- 
ers with more than two progenitors, or mass accretion 
including all parents. However, when considering short 
timescales (< 1 G yr), the two body c riterion is a good 
approximation (see lWetzel et al.l (|2007f ) for more details). 

All errors cited are la errors derived from dividing 
the simulations into 8 octants and computing the dis- 
persion of the quantity of interest within each octant. 
Since we probe scales much smaller than the octants, we 
treat them as uncorrelated volumes. 

3. CLOSE PAIRS AS PREDICTORS OF MERGERS 

To examine whether close pairs of galaxy clusters form 
reliable predictors o f mergers, we first extend the work of 
IWetzel et ail (|2007f ) on halo mergers, using the same sim- 
ulation and halo catalogs (described in $2]). We identify 
child halos of mass M > 5 x 10 13 }i~ 1 Mq that are prod- 
ucts of recent mergers with progenitor masses M\ and 
M2, where M2/M1 > 0.2 (major mergers). We explore 
the distribution of progenitor separations in Fig.[TJ which 
shows the number density of progenitors as a function 
of binned progenitor separation. For time intervals of 
1 Gyr, there is a characteristic comoving progenitor sep- 
aration (~ 1.5 /i _1 Mpc) that does not evolve significantly 
with redshift. Without the influence of gravitational at- 



traction, halos with typical velocities of ~ 1000 km/s will 
travel 1 Mpc within 1 Gyr. Furthermore, the number 
density of mergers evolves only weakly with redshift. 
These factors suggest that close pairs might be a rea- 
sonable proxy for mergers. 

These results represent a post-diction, where we know 
in advance which halos will merge. We now investigate 
whether close pairs of objects can reliably be used to pre- 
dict mergers. The following subsections investigate how 
frequently pairs will merge as a function of pair mass, (co- 
moving) separation, redshift, and merger time interval. 
We consider a given pair to have merged if both halos are 
progenitors of the same child at a later timestep. Since 
observations of galaxy clusters are generally limited by a 
threshold luminosity, we identify pairs of halos whose in- 
dividual halo masses are above a given mass cut. We first 
examine the pair-merger hypothesis (that close pairs are 
merger proxies) for the best-case scenario in which halo 
positions and masses are known with complete accuracy 
(as in our simulations). This will provide a firm upper 
limit to the utility of pairs in predicting mergers. Then in 
2] we consider observational complexities such as scatter 
in the halo mass, redshift space distortions, and redshift 
space errors, which significantly degrade the signal. 

3.1. Pair Mergers at z = 

We select halos at z — 0.04 and consider mergers to 
z = (At = 0.6 Gyr). Figure (top) shows the number 
density of pairs as a function of pair separation, n p (r), 
along with the number densities of pairs that merged 
(solid), n m (r), and did not merge (dashed) within the 
time interval. The number density of pairs terminates at 
small separations because of halo exclusion, i.e. massive 
halos have finite radii. The upper set of curves are halos 
of mass M > 5 x 10 13 h^ 1 M Q and the lower have M > 
10 14 h~ 1 M & ; the results are qualitatively similar. There 
is a limited range of separations larger than the halo 
exclusion limit in which the majority of pairs merge. Also 
plotted (center) are the same number densities, but as 
an integrated function of separation, demonstrating the 
relationship for pairs within the given separation r. In 
0.6 Gyr at z = 0, massive halo pairs will only merge with 
certainty for separations < 2 /i~ 1 Mpc. 

Figure [2] (bottom) also shows the fraction of pairs that 
merge within the given separation, n m {< r) / n p {< r) 
(falling curves). This fraction is the "efficiency" of the 
close-pair method in identifying merger candidates, since 
it shows the likelihood that a pair within the given 
separation will merge. It can also be thought of as a 
measure of contamination of the candidate sample, be- 
cause its difference from 1 identifies the fraction pairs 
that do not merge. The rising curve shows the frac- 
tion of mergers over the total number density of mergers, 
n m {< r)/n m j t- This can be interpreted as the merger 
"completeness," as it shows the fraction of the total num- 
ber of mergers found by identifying pairs within the given 
separation. Figure [2] confirms the intuitive result that 
because of greater accelerations, larger halos are able to 
merge from larger separations, which holds at all red- 
shifts and time intervals (see Fig. [3}. 

If close pairs were ideal predictors of mergers, the ris- 
ing and falling curves of Fig. [2] (bottom) would be step 
functions, crossing near a fraction (y-value) of 1 and 
demonstrating a clear dichotomy between those pairs 
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M cut =5e13 (upper), 1e14 (lower) Mg/h 
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Pair Separation r (Mpc/h) 



Fig. 2. — Top: The number density of pairs at z = 0.04 as 
a function of binned separation, n p (r) (black solid), the subset of 
these that merged in At = 0.6 Gyr, n m (r) (blue solid), the remain- 
der that did not merge (red dashed). The upper curves (thin) are 
for pairs of mass M > 5 X 10 13 h~ 1 Mq, and the lower (thick) are 
for more massive pairs, M > 10 14 /i _1 Mq. The intersection with 
the dashed curve shows the separation beyond which the majority 
of pairs do not merge, highlighting the limited range of separations 
for which merging is likely. Center: The same, but integrated for 
all separations smaller than the abscissa. Bottom: The fraction 
of pairs separated by less than r that merged, n m (< r)/np(< r) 
(falling curves), indicating the efficiency of pairs as merger can- 
didates, and the ratio of identified mergers to all mergers above 
the mass cut, n m (< r)/n m ,tot (rising curves), indicating the com- 
pleteness of the candidate sample. Error bars and shaded regions 
indicate the lcr error of cosmic variance from dividing the simula- 
tion into octants. 

that merge and those that do not. However, because 
of dynamical effects halos do not simply accelerate to- 
ward one another from rest; they merge with neighbors 
from a broad distribution of separations. Indeed, the ef- 
fects of local environment can hinder the merging of close 
pairs that would ordinarily occur via non-linear two-body 
gravitational interaction (see [j5]for more details). 

The intersection of the efficiency and completeness 
curves represents an easily identifiable pair separa- 
tion that compromises between maximizing completeness 
while minimizing contamination. At this intersection the 
fraction of pairs that merge, n rn (< r)/n p (< r), equals 
the fraction of total mergers, n m {< r)/n m ^ot, although 
the relevant populations are different. If this intersec- 
tion occurs at a fraction close to 1, most pairs within the 
separation will merge, and those pairs will represent the 
majority of all mergers that occur in the time interval. 
However, if this fraction is low, pairs within the separa- 
tion are not a representative indicator of mergers. Since 
the intersection of these two curves occurs at a fraction 
of ~ 0.5 for massive pair mergers within At = 0.6 Gyr 




Physical Pair Separation r (Mpc) 

Fig. 3. — The fraction of pairs separated by less than r that 
merged (falling curves), indicating efficiency, and the ratio of iden- 
tified mergers to all mergers above the mass cut within At = 
1 Gyr (rising curves), indicating completeness, for various red- 
shifts. The total number density of such mergers is 458 h~ 3 Gpc -3 
(z = 0.97-0.75), 950/i~ 3 Gpc- 3 (z = 0.58-0.44), 1464 h" 3 Gpc- 3 
(z = 0.12 — 0.04). Note the approximate invariance of the total 
merger fraction (rising curve) with redshift. Bottom: The same, 
but as a function of physical separation. While pairs of a given 
comoving separation merge more frequently at high redshift, this 
trend reverses when considering physical separations. 

at z = 0, we conclude that pairs at these redshifts and 
masses cannot be used to reliably predict mergers. 

3.2. Redshift Dependence of Pair Mergers 

We next examine whether the pair-merger hypothe- 
sis fares better at high redshift. We consider a longer 
merger timescale of At = 1 Gyr for massive halo pairs 
at z a 1,0.6, and 0.1. Figure [3] (top) shows the pair 
merger fractions as a function of (comoving) separations. 
The pair merger efficiency (falling curves) increases sig- 
nificantly with redshift, while the merger completeness 
(rising curve) remains approximately invariant. For a 
given merger timescale, mergers come from pairs of es- 
sentially the same (comoving) separation, regardless of 
redshift (see Fig. [T]). However, since gravitational inter- 
actions are governed by physical (not comoving) separa- 
tions, Fig. [3] (bottom) shows the same merger fractions 
as a function of physical separation. The merger effi- 
ciency trend is reversed: pairs within a given physical 
separation are more likely to merge at low redshift, and 
mergers come from pairs of much larger physical separa- 
tions at low redshift. The intersection of the efficiency 
and completeness curves occurs at a higher fraction at 
high redshift than at low redshift (75% vs 60%), indicat- 
ing that, for a fixed merger timescale, massive halo pairs 
provide a slightly better proxy for merger rates at higher 
redshift. 
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Fig. 4. — Top: The fraction of pairs separated by less than 
r / r vir that merged (falling curves), indicating efficiency, and the 
ratio of identified mergers to all mergers above the mass cut within 
At = 1 Gyr (rising curves), indicating completeness, at z ss 1 
for various mass cuts. The total number density of such merg- 
ers is 45,750fe- 3 Gpc- 3 (M > 5 x 10 12 /i- 1 M©), 3957 /i^Gpc -3 
(M > 2 x 10 13 h~ 1 M Q ), 458h" 3 Gpc- 3 (M > 5 x 10 13 ft- 1 M ). 
Middle: The same, but for for higher mass halo pairs merg- 
ing within At = 0.5 Gyr of z = in the (2.2 h _1 Gpc) 3 simula- 
tion. The total number density of such mergers is 1000 /i~ 3 Gpc — 3 
(M > 5 x IO^/i^Mq), 280h" 3 Gpc- 3 (M > 1 x 10 14 h~ l M Q ), 
56/i -3 Gpc- 3 (M > 2 x 10 14 h- 1 M Q ). Bottom: The same as 
middle, but for a linking length of b = 0.20. Masses and virial 
radii have been scaled to that of b = 0.15 for direct comparison. 

3.3. Mass and Linking Length Dependence of Pair 
Mergers 

While it appears that cluster-mass halo pairs might 
provide a reasonable proxy for mergers at z ~ 1, we next 
examine whether these results are robust as a function 
of halo mass. Figured] (top) shows the pair merger frac- 
tions for halo masses down to M > 5 x 10 12 h~ 1 M Q at 
z Ri 1, where the pair separation is shown as a multiple 
of the minimum mass cut virial radius (r^cio). Figured 
(middle) shows the same but for higher halo masses and 



a shorter time interval at z w 0. While pairs of higher 
mass halos merge from larger physical pair separations, 
when scaled by the halo virial radius we find that the 
pair merger efficiency and completeness exhibit a nearly 
universal relation as a function of halo mass. This im- 
plies that any of our results in for a given mass can 
be approximately scaled to those of another mass. 

This universal merger relation also implies that the to- 
tal merger fraction can be estimated by noticing that 
the following relation is satisfied at the crossing point 
of the curves: (n m (< r)/n p {< r))/(n m {< r)/n„ htot ) = 
n m ,tot/n p {< r) = 1. For example, from the middle panel 
of Fig. U we find that ^ m) tot = n p{< r ) when evalu- 
ated at a scaled pair separation of 5r v j r . One can thus 
estimate the merger rate by counting pairs interior to 
5r v j r . This is a potentially powerful result since it is 
approximately invariant for all mass-scales. 4 However, 
this radius is a function of the timescale considered, and 
is likely to depend on the cosmological parameters. Fur- 
thermore, although their number densities agree, the pair 
population within 5r vn , does not directly corespond to 
the total merger population. This makes it impossible to 
use the pair population to directly study any other prop- 
erties of mergers (e.g., clustering). Redshift space dis- 
tortions will significantly undermine the utility of these 
results (see©. 

The only strong deviation from this nearly univer- 
sal mass relation occurs for merger completeness of 
5 x 10 12 /i _1 M Q halos (top), which results in an inter- 
section of merger efficiency and completeness at a lower 
fraction (10-15% decrease) for lower-mass halos, a trend 
which we find does not depend strongly on redshift or 
merger timescale. The use of pairs as proxies for merger 
rates becomes increasingly unreliable as we approach 
massive galaxy-size halos, an important result for mas- 
sive galaxy mergers at high- redshift, where galaxies are 
found primarily in distinct host ha los. 

These results augment those of iBerrier et al.l (|2006f l. 
who found that the evolution of close galaxy pairs cannot 
be used to measure the host halo merger rate. Specifi- 
cally, they found that the host halo merger rate evolves 
as (1 + z) 3 , while the number of close galaxy pairs evolves 
little with redshift. This arises because, at low redshift, 
the merger rate of host halos is low, but there are multi- 
ple galaxy pairs merging within massive host halos. At 
high redshift, the host halo merger rate is high, but the 
number density of halos massive enough to host more 
than one bright galaxy is low. Our results imply that 
even when considering major mergers of massive galax- 
ies at high redshift, where each halo hosts a single bright 
galaxy, the close pair population will not reliably trace 
the merger population. 

We have also compared two different linking lengths, 
to explore the dependence of our results on the FoF pro- 
cess. In Fig. |4] (bottom), we show the results using a 
linking length of b — 0.20. Since changing the linking 
length changes the mass of a given halo, we have scaled 
the mass threshold for b = 0.20 to match the number 
densities of the b — 0.15 sample, thereby probing the 
same halo population. We have also scaled the halo sep- 
arations by r v j r = rsoo- In the case of b = 0.20, the 

4 We thank the referee for pointing out this universal behavior. 
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Pair Separation r (Mpc/h) 



Fig. 5. — The fraction of pairs separated by less than r that 
merged (falling curves), indicating efficiency, and the ratio of iden- 
tified mergers to all mergers above the mass cut (rising curves), in- 
dicating completeness, for various merger timescales from z = 0.58. 
The total number density of such mergers is 950/i" 3 Gpc- 3 (z = 
0.58 - 0.44), 2755fe- 3 Gpc- 3 (z = 0.58 - 0.22), 4282 /i^Gpc" 3 
(z = 0.58 - 0.04). 

intersection of the efficiency and completeness curves is 
shifted outward by ~ r v ^ T , which is not surprising since 
a halo found using b = 0.20 is expected to be ~ 30% 
larger than one found using b — 0.15. From the figure we 
see that changing the linking length results in a change 
of only a few percent in the merger fraction at the inter- 
section of the efficiency and completeness curves. Thus 
changing the linking length does not qualitatively alter 
any of our results. The weak dependence on the linking 
length suggests that for these mergers, the FoF proce- 
dure does not give rise to significant artificial bridging 
of nearby halos, and our results would not differ signifi- 
cantly from a similar an alysis performed usin g a spherical 
overdensity halo finder (|Lacev k, Cold [19941 ). 

3.4. Merger Timescale Dependence of Pair Mergers 

Pair merging is also strongly dependent on the choice 
of merger timescale. Figure [5] shows that the merger effi- 
ciency (falling curves) increases with timescale, i.e. pairs 
are more likely to merge when one considers longer 
time intervals. When considering mergers across 5 Gyr, 
a significant fraction (80%) of pairs with separations 
< 5/i _1 Mpc will merge. However, the merger complete- 
ness within a given separation (rising curves) decreases 
with timescale, so one must consider pairs at larger sep- 
aration to capture all mergers. Interestingly, the trade 
off between efficiency and completeness is only weakly 
dependent on the timescale, and we find similar results 
even for massive halo pairs at z = 1 merging to z = 
(At = 7.6 Gyrs). No matter which merger timescale 
are considered, the utility of pairs as proxies for mergers 
remains limited. 

In all redshift, mass, and temporal regimes we con- 
sider, close halo pairs do not provide a robust predictor 
of overall merger rates. We conclude that, while pairs 
at small separations can reliably predict mergers, these 
constitute a small fraction of the total merger popula- 
tion. While our results could in principle be used to cali- 
brate the merger fraction as a function of halo separation, 
close pairs do not provide a self- consistent probe of the 
merger population independent of our theoretical predic- 
tions (which are contingent upon having simluated the 



correct cosmology). Thus, even in the best-case scenario 
of complete knowledge of mass and position, measure- 
ments of galaxy cluster pairs cannot be used to measure 
cluster merger rates. The same is true for any any class 
of object that singly occupies a halo, particularly in using 
field galaxy pairs at high redshift to probe galaxy mergers. 

4. SCATTER IN MASS, REDSHIFT SPACE DISTORTIONS, 
AND REDSHIFT ERRORS 

The situation becomes more complicated when we con- 
sider redshift space distortions, redshift space errors, and 
scatter in the estimated mass of the halos. To identify 
the impact of imprecise determination of dark matter 
halo masses, we have recomputed several of the statistics 
of after introducing an RMS scatter of 0.2A/ cut , where 
M cut is the threshold for detecting the pairs in the mock 
observation. The scatter in the mass causes some halos to 
fall out of the sample and others to enter it, resulting in a 
different population of halos near the threshold. Because 
the mass function is steep, many more low-mass halos en- 
ter the sample than high-mass halos leave it. Thus there 
are both more pairs to consider as merger candidates and 
more actual mergers between members of the observed 
sample. For example, if M cut = 5 x 10 13 ft _1 M Q , the 
number of pairs increases more than the number of merg- 
ers, resulting in a decrease in n m /n p of ~ 5% for pairs 
separated by 2 < r < 6/i _1 Mpc over 0.58 < z < 0.97 
(At = 2 Gyrs). The effect is similar for other thresh- 
olds, redshifts, and intervals. Scatter in the mass has a 
more pronounced impact at high values of M cut where 
the mass function is steeper. The trend with redshift 
is derivative of the evolution of the mass function: at 
a fixed value of M cut there is more sensitivity to scat- 
ter at higher redshifts where the mass function steepens. 
There is not a significant trend with merger timescale 
At. The overall impact of imprecise mass measurements 
is to degrade the range in separation over which pairs 
might be considered useful proxies for mergers, but it is 
not a significant effect. 

In contrast, redshift space distortions from doppler 
shift due to the halo peculiar velocities v p are catas- 
trophic for the pair-merger hypothesis. When pairs are 
identified in redshift space, virtually none of them merge. 
In redshift space, the line of sight component is shifted 
relative to its real position by an amount 

r*+*,/c cdz 

Al " - 1 m (1) 

For typical halo velocities in our simulation, this shift 
amounts to a few h~ 1 Mpc. Figure [6] illustrates that this 
distortion makes many close pairs appear to be highly 
separated and pairs that are highly separated to appear 
close. At redshift space separations less than 5/i _1 Mpc, 
fewer than 5% of merger candidates actually merge in 
At = 0.6 Gyr at z = 0. Even at z = 0.58 where redshift 
space distortions are less severe, only ~10% of candidates 
merge over At = I Gyr. The situation improves modestly 
for longer time intervals, with ~60% of pairs merging 
over At = 5.6Gyr, but this falls far short of the fraction 
in real space. 

Figure [5] shows a histogram of the difference between 
a pair's separation in redshift space (dashed line), as 
would be done in an observation, and its true separa- 
tion in configuration space (solid lines) at z = 0.04. 
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Fig. 6. — A histogram of the change in apparent separation 
because of redshift space distortions from halo peculiar velocities. 
Dashed lines indicate close pairs identified in redshift space (as in 
observation) and solid lines indicate pairs identified in real space. 
The region < of the dashed lines indicate pairs that are artificially 
close in redshift space. The region > of the solid lines indicate 
pairs that are artificially separated in redshift space. 

Merger candidates identified in real and redshift space 
are widely disjoint sets. The three sets of curves (colors) 
are bins of separation, as annotated. For the population 
of merger candidates identified in redshift space (dashed 
lines), the region greater than zero shows pairs for which 
the real separation was smaller than the measured sep- 
aration, while the negative region shows pairs that ap- 
peared to be close in redshift space but were actually 
separated by several ft,~ 1 Mpc. The former set might be 
likely to merge, but the latter set is highly unlikely to do 
so. For the actual population of close pairs (solid lines), 
the region greater than zero shows pairs whose redshift 
space separations were greater than their real separa- 
tions, and which were therefore likely removed from the 
sample of merger candidates. The reason for the asym- 
metry between the pair population that left the sample 
(solid, right side) and the population that entered the 
sample and replaced it (dashed, left side) is that pairs 
are most often found in overdense regions, toward which 
many other halos are streaming. The impact is that the 
pair population identified in redshift space is larger and 
has a broad distribution of physical separations; virtu- 
ally none of the "merger candidates" we identified actu- 
ally merge. This result is robust for all redshifts, merger 
intervals, and halo masses. 

To a degree, this result also applies to the relationship 
between galaxy pairs and galaxy mergers. Observing 
galaxy pairs to study the merger rates of singly occu- 
pied host halos at high redshift will be similarly affected 
by redshift space distortions, as will using galaxy pairs to 
deduce the merger rates of multiply occupied host halos 
at lower redshift. However, the impact of redshift space 
distortions on using galaxy pairs to deduce the galaxy- 
galaxy merger rates within a single host halo may be less 
severe (jBerrier et al.l 1200(f ). Inside a halo, galaxy merg- 
ers are predominately between a satellite galaxy and the 
central galaxy, and merger dynamics inside a host halo 
are influenced by dynamical friction, especially in group- 
mass halos and when the galaxy masses are comparable. 
Pair counts, even in redshift space, provide a rough proxy 
to the halo mass and will therefore correlate well with the 
merger rate within the halo. 



The impact of redshift space errors is similar to that 
of redshift space distortions, except that, unlike with the 
velocities, the errors are not correlated with the high 
density regions. Redshift space errors that result in more 
than a few /i _1 Mpc shift in the apparent position have a 
similar impact to redshift space distortions. 

5. THE MERGER KERNEL AND THE DENSITY 
DEPENDENCE OF MERGERS 

Since spatial information is a weak probe of halo 
merger statistics, we now turn to statistically describing 
the merger rate via the halo mass function. For any 
population of objects that are built up by binary mergers 
of smaller constituents, the rate at which the number 
of objects of a given mass changes can be described 
by the Smoluchowski coagulation equation, which has 
been applied to th e evolution of dark matter halos 
(ISilk fc White) Il978t ICavaliere. Colafrancesco. fc Mencil 
1992 ~lSheth fc Pitmanl 119971 : 

Benson. Kamionkowski fc Hassanil 120051 ) . The abun- 



dance of halos at a given mass is increased through the 
creation of such halos by mergers of smaller halos, and 
decreased as these halos merge into more massive ones. 
The rate of change of the number density of halos of 
a given mass can then be determined by knowing the 
number density of halos at all masses, i.e. the mass 
function, and the proper merger kernel to relate the 
mass function to a merger rate. Historically, it has 
been difficult to study the coagulation of cluster-mass 
halos through simulation because these events are rare. 
However, with our large simulation volume the merger 
kernel can be computed in a statistically significant way. 

We define the merger kernel (J(mi,m2, z) as in 
iFurlanetto fc Kamionkowski (|2005l ) (hereafter F&K), via 
the relation 



Q(mi,m 2 , z) 



n m (m 1 ,m 2 , z) 
n(mi, z) n(rri2, z) At' 



(2) 



where n m (m\, vn 2 , z) is the number density of merg- 
ers between parents of mass mi and m 2 in a time At, 
and n(m, z) is the halo mass function. The quantity 
Q can be interpreted as the efficiency of merging for a 
pair of objects, such that the rate of mergers is a prod- 
uct of this efficiency with the densities of available par- 
ents. Note that in previous contexts the term efficiency 
refers to the ability of the close-pair technique to re- 
liably identify merger candidates, whereas in this con- 
text the efficiency is the coefficient of the number den- 
sities when quantifying the merger rate. The time in- 
terval At should be sufficiently short that there is no 
significant evolution of the halo mass function. Satisfy- 
ing this requirement drives down the number of mergers 
in a fixed volume, making statistically significant mea- 
surements difficult. However, with our large simulation 
volume we measure Q(m\,m2, z — 0) on the interval 
< z < 0.04 (At = 0.6 Gyr) by classifying the two pro- 
genitors that contributed the most mass to each halo at 
z = into ten logarithmically spaced mass bins in the 
range 10 13 < h~ 1 M Q < 10 15 . The results are shown in 
Fig.0 

We find that Q follows the simple functional form 



Q(mi,m 2 , z = 0) = A 



mi + m 2 



1 B 







(3) 
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TABLE 1 

Best fit amplitude, A, and exponent, B, for the merger 

KERNEL Q = A(m 1 + m2) B , WHERE (mi + 1712) IS THE SUM OF 

THE 2 PROGENITOR MASSES. THE FIRST FOUR ROWS ARE THE 
RESULTS FOR EACH OF THE 4 DENSITY SUBDIVISIONS, WHILE THE 
BOTTOM ROW IS THE BEST FIT FOR THE ENTIRE SIMULATION. THE 
AVERAGE DENSITY AND TOTAL VOLUME OF EACH SUBDIVISION IS 
LISTED AT LEFT, AND THE REDUCED \ 2 F0R EACH SUBDIVISION IS 
LISTED AT RIGHT. 
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The best fit values for A and B, found using a lin- 
ear least-squares fit to the data, are presented in 
Table Q] (bottom row). This functional form satis- 
fies the formal requirement that the merger kernel 
be s ymmetric in its two arguments. For compari- 
son, iBenson. Kamionkowski fc Hassanl (|2005f ) analyti- 
cally found that Q oc (mi + m%) when P[k) oc fc™ with 
n = 0, which is approximately true in the trans- linear 
regime of cluster-mass halos. 

5.1. Density Dependence of the Merger Kernel 

Because the densities of the parent populations have 
been divided out, it is conceivable that the merger kernel 
Q is independent of the large-scale density field. Signifi- 
cant dependence on density would indicate that environ- 
mental effects other than the progenitor densities, such 
as halo velocity distributions, halo impact parameters, 
and tidal fields, are at play in determining the merger 
rate. To investigate the dependence of Q on the large- 
scale density field we construct a coarse (64 3 ) density 
grid in our simulation and assign the dark matter parti- 
cles to the nearest grid point, effectively smoothing the 
field on a scale of ~ 17ft, _1 Mpc. Each halo at z = 
is also assigned to the coarse grid, and these halos are 
sorted in order of their large scale density environment 
and then divided into quartiles of density, with each quar- 
tile containing approximately the same number of halos. 
Using the coarse grid, we compute the total volume and 
overdensity of each quartile in the simulation, and these 
appear in the first two columns of Table [TJ The mean 
overdensity of a quartile is defined as 

(Si) = {p)i ~ {P) i = 1,2,3,4 (4) 

(P) 

where the angular brackets denote an average of the 
~(17 /i _1 Mpc) 3 cells of the coarse grid. The merger ker- 
nel Q is fit for each quartile to the form of Eq.02 using the 
value of B determined from the entire simulation. The 
results are summarized in Table [T] The best fit values of 
A differ by several sigma, and the improvement in the re- 
duced Xicd statistic indicates that each of the individual 
density fits is a much better fit that a fit to the entire vol- 
ume. This indicates that there is a clear trend in Q with 
the large-scale density field: as the density increases the 
efficiency of merging for a given system mass decreases. 
Thus, the environmental effects in dense environments 
are hindering the merging process in comparison to less 
dense environments. 



Fig. 7. — Merger kernel Q from the simulation for 10 logarithmic 
mass bins in progenitor masses, as a function of the sum of the 
progenitor masses. The solid line depicts the best fit of Eq. [3] to 
the whole simulation volume. The shaded region corresponds to 
the formal la errors on the fit values of A and B. The top panel 
shows Q computed in the highest and lowest density quartiles and 
the bottom panel shows the quartiles closer to the mean density. 
The highest and lowest density regions differ in Q by several sigma, 
and the density dependence of Q is one of the main factors driving 
the dispersion in the fit to the whole box (shaded region). 

Figure [7] shows Q in each of the four density environ- 
ments plotted individually. The upper panel shows Q in 
the highest and lowest density quartiles, while the lower 
panel shows the results from the inner two quartiles. The 
solid line is the best fit model for the entire box, which 
is only a good approximation for environments close to 
the mean density. The yellow shaded region results from 
allowing A and B to deviate by their la values. Figure[7] 
demonstrates clearly that a large component of the dis- 
persion in the best fit to the entire simulation originates 
in the density dependence of Q. 

Figure [7] clearly indicates that merging is less efficient 
in high density environments than in low density envi- 
ronments, but we note that with regard to our fit there 
is an ambiguity as to whether this implies a lower nor- 
malization A or a shallower power law B. Given the 
rarity of these events at the high mass end and the lim- 
ited dynamic range that is consequently driving the fit, 
there is large covariance between the fit values of A and 
B. We have performed the 2-parameter fit to Eq. [3] for 
each of the density quartiles individually and find that in- 
deed A increases slightly with density (in contrast to the 
1-parameter fit) while B decreases. Unfortunately, the 
formal errors in the parameters (from inversion of the 
covariance matrix) become so large that the four regions 
are statistically indistinguishable. However, the individ- 
ual fits hint at an interesting possibility: if B decreases 
and A increases with density, then the merger efficiencies 
cross over as a function of system mass. Specifically, for 
lower total system masses (mi + TO2), merging becomes 
more efficient in high density environments. In our sim- 
ulation, all four Q curves cross over in the mass range 
\Q w h~ l M s < (mi + ma) < 10 12 /i _1 M Q , suggesting 
that galaxy scale mergers may be more efficient in denser 
environments. While this evidence is extrapolated from 
higher mass, the trend is potentially worth further inves- 
tigation given that other halo properties, e.g. clustering 
as a function of concentration/ formation time, reverse 
trend from M > M* to M < Af* (jGao, Springel fc White! 
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Fig. 8. — The fraction of the total number of pairs sepa- 
rated by less than r that merged (falling curves), indicating ef- 
ficiency, and the fraction of all mergers in the simulation within 
At = 1 Gyr (rising curves), indicating completeness, for the high- 
est and lowest density quartiles. Pairs are found at z = 0.12 for 
M > 5 X 10 12 h~ 1 M e (top) and M > 10 14 h- 1 M @ (bottom). 



[20051 IWechsler et alJl2006l) 5 . 

We note that the merger kernel Q as computed in this 
section is not a direct cosmological observable. There is 
an observational counterpart to the merger kernel: the 
pair kernel, Q p , computed using the number density of 
pairs with separations r < r& t - The threshold separation 
is calibrated from simulations to any desired tolerance for 
completeness or contamination, but the result is insensi- 
tive to r-At- In all cases, there are many more pairs than 
actual mergers, the amplitude A of Q p is larger, and the 
power law B is shallower, both by several a. It infeasi- 
ble to use pairs at any mass scale as a proxy to compute 
merger rates. 

5.2. Density Dependence of Close Pair Mergers 

The trend in the merger kernel is driven by the fact 
that the number of mergers grows more slowly with 
density than the parent mass functions. By studying 
changes in the relationship between close pairs and merg- 
ers in regions of differing local density, we can demon- 
strate that this is not universally true for all pair sepa- 
rations. We proceed in an analogous manner, smooth- 
ing over ~ 17/i _1 Mpc and defining four density quar- 
tiles, each containing the same number of mergers. We 
consider pairs merging between z — 0.12 and z — 0.04 
(At = 1 Gyr), though our results are insensitive to 
the time interval and redshift. Fig. [5] shows the re- 

5 M„(z) is the mass at which cr(M), the variance of the linear 
power spectrum smoothed on scale M, equals the threshold for 
linear density collapse 8 c (z); A/* ~ 8 X 10 /i _1 Mq at z = 0. 



suits for pairs of halos above 5 x 10 12 h~ 1 M Q (top) and 
10 14 h~ 1 M & (bottom). The relationship between close 
pairs and mergers changes with the large-scale density. 
Pairs within a given separation are a less complete sample 
of merger candidates in overdense environments than in 
underdense regions (rising curves, top and bottom). Ob- 
jects in overdense regions have higher velocities, which 
allow mergers to come from a broader range of progen- 
itor separations and can extend the infall time of closer 
pairs by generating larger angular momenta. 

Figure[S]also indicates that, at large separations, pairs 
are more efficient predictors of mergers in overdense re- 
gions than in underdense regions (falling curves) because 
higher velocities in overdense regions allow pairs from 
larger separations to merge. However, at small separa- 
tions (< 2/i~ 1 Mpc), the relationship between pairs and 
mergers depends on the pair masses. Here, there are no 
high-mass (> 10 14 /i~ 1 Af Q ) halos because of halo exclu- 
sion, and pairs are more efficient predictors of mergers 
in underdense regions (top). This is likely because low- 
mass progenitors in overdense regions have large veloc- 
ities and thus have longer infall times, e.g. from tidal 
fields, resulting in large angular momenta that inhibit 
merging. In contrast, the smaller velocities in underdense 
regions enhance the probability that a very close pair 
will merge from simple gravitational attraction. High- 
mass halos have much stronger gravitational attractions 
and are thus less affected by the large-scale density field. 
Indeed, at small separations, massive halos are equally 
likely to merge in over- and underdense environments. 

Finally, we note that the two methods we have pre- 
sented to examine the density dependence of merging 
probe mergers in different ways. The merger kernel in- 
tegrates over merger pair separations, thus considering 
mergers as a function of progenitor pair mass sum. Al- 
ternately, the close pair merger method selects only close 
halo pairs that are above a given mass cut, thus integrat- 
ing over mass dependence and considering mergers as a 
function of pair separation. The merger kernel indicates 
that for all child halo masses in our simulation merg- 
ing is more efficient in underdense environments. This 
is corroborated by Fig. [8] (top), which shows that for 
all halos in our simulation, pair merging is more likely 
in underdense environments for pairs within 1.5 h~ 1 Mpc 
separation. The accompanying merger completeness (ris- 
ing) curves show that pairs within this separation consti- 
tute the majority of pair mergers in the given timescale. 
Figure [8] (bottom) shows that high-mass pairs are more 
likely to merge in overdense regions at all pair separa- 
tions. However, since mergers between two halos both 
with M > 10 14 h~ 1 M ( r ) are extremely rare (as opposed 
to mergers between a high and low mass object), these 
contribute less to the merger kernel at a given child mass. 



6. 



THE CLUSTERING OF CLOSE PAIRS AND MERGER 
BIAS 



For two halos to merge, they must have been located in 
close physical proximity at an earlier time. Since closely 
spaced halos are more likely to be found in overdense re- 
gions, recently merged halos may exhibit enhanced spa- 
tial clustering. Moreover, if recently merged halos cluster 
differently from the general population ( "merger bias" ) 
and this is unaccounted for, conclusions drawn about 
halos on the basis of their clustering could be compro- 
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Wetzel ct al. 



mised. For example, the use of cluster self-calibration 
to infer cluster masses (and in turn cosmological pa- 
rameters) requires a precise knowl edge of the cluster- 
ing o f clusters as a function of mass (jMaiumdar fe Mohrl 
2003). A number of earlier studies have looked for such 
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mixed result^ 

Recently, iFurlanetto fe KamTo nkowskij (|2005l ) devel- 
oped an analytic model to predict the merger bias. As- 
suming that mergers correspond to closely spaced objects 
at an earlier time, they compare the clustering of close 
pairs to that of single objects, thus computing the pair 
bias as a proxy for merger bias. On scales much larger 
than the pair separations, they found an enhancement 
of clustering for pairs of mass M > M*, implying that 
recently merged high-mass halos should exhibit a clus- 
tering bias. We now use F&K's framework to determine 
if the clustering of close halo pairs of mass M ^> M„ 
provides an accurate p roxy for the cl u sterin g of mergers. 

Using simulations, IWetzel et al.l (|2007f) found the 
most prominent merger bias (~ 20%) for major merg- 
ers (M 2 /Mi > 0.3) of high-mass halos (M > 5 x 
lO 13 /i _1 Af0) at z = over At = 0.6 Gyr. We con- 
sider a similar mass and temporal regime. To improve 
our statistics we additionally use a box of eight times the 
volume previously used (see S}2]) , which allows us to probe 
merger effects of more massive halos. We use a shorter 
timescale of At = 0.5 Gyr at z = to define our merger 
interval to preserve the signal in the comparison popula- 
tion. Contrary to expectation, the merger bias does not 
increase with mass; it remains a 10-20% enhancement 
with similar statistical significance up to halos of mass 
M > 4 x 10 14 h~ l MQ. 

In computing the pair bias, F&K define "close" halo 
pairs by demanding that the probability of finding three 
or more halos in a sphere of a given radius is small 
compared to that of finding two. This is approximately 
equivalent to the restriction that the probability of find- 
ing two or more neighbors, P(~> 2), within a given sep- 
aration from a halo is small compared to that of finding 
one, -P(l), which is what we will measure. For halos M > 
5xl0 13 h~ l M Q atz » 0, we find that P(> 2)/P(l) » 0.1 
for a (comoving) sphere radius of 4 /i -1 Mpc. This separa- 
tion restriction yields ~ 6,800 pairs per (ft, _1 Gpc) , out 
of ~ 77,000 halos per (/i _1 Gpc) 3 . While this is a suf- 
ficient number density of pairs for a robust correlation 
function measurement in our (2.2 ft, Gpc) simulation 
volume, [^demonstrated that close pairs do not reliably 
predict the merger population. 

The analytical model of F&K predicts a significant 
merger bias in its application to the clustering of mas- 
sive halos (M ^> M*). For such objects, it predicts a 
correlation function of pairs X p (r) in terms of Oj(r), the 
correlation function of individual halos: 



X p (r) = [l + Ur)} 4 



1. 



(5) 



This is computed on scales where the underlying matter 
fluctuations are linear, which at z = corresponds to 
distances greater than ~ 10/i _1 Mpc. To compute the 
pair bias F&K use 
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Fig. 9. — Correlation function of couples (triangles) at z = 0.04 
compared with that of single halos (squares) at z = 0. Mass cuts 
are chosen so that the single halos are approximately the same mass 
as the child products of the couples. While enhanced signal is seen 
at 47h~ 1 Mpc, this is mitigated by the adjacent points which are 
both negative. 



For a given halo mass, this can result in anomalously 
high pair biases. This is because the halo correlation 
function is implicitly a function of the halo mass; more 
massive halos are more highly clustered. When com- 
puting the merger bias (or its proxy, the pair bias) the 
comparison population should be an ensemble of halos of 
the same characteristic mass as the child halo, not the 
parent halos. Otherwise, the merger bias becomes entan- 
gled with the effect that larger halos are more clustered. 
To adjust for this in the analytic model we note that 
&(r,2M) « l.l&(r,M) near M ~ 10 14 M Q . On scales 
> 10/i _1 Mpc, £(r) is less than 1, so this leads to a pair 
bias of bp » 3.6. This pair bias is still significantly larger 
than the 10-20% merger bias seen in simulation for halos 
up to M >4 x lQ 1A h- l M Q . 

As an alternative to an analytic approach, we next 
measure the clustering of close pairs in simulations to 
discover whether pair bias can predict merger bias. We 
consider pairs 0.5 Gyr prior to z = 0. To assign a 
unique position to a pair of neighboring halos, we im- 
pose a "couples" restriction, namely that each member 
of a pair is the other's nearest neighbor. This restriction 
remains robust for two body mergers, a good approxima- 
tion for our time interval. Within the pair separations we 
consider (< 2 h~ 1 Mpc) couples constitute essentially all 
pairs. We select couples of halos above 5 x 10 13 h^M® 
at z = 0.04 and use their geometric centers to evaluate 
the couple correlation function. We use couples with sep- 
arations less than 1.6/i _1 Mpc, 80% of which correspond 
to mergers in our time interval (Fig. [J). To limit the ef- 
fect of mass-dependence on the halo correlation function, 
we compare this correlation function with that obtained 
from sing le halos above 10 14 /i^M© at z = 0. While 
this ignores the effects of mass scatter in mergers, such 
effects remain small for this short timescale. 

Figure [9] shows that with a (2.2 /i _1 Gpc) 3 simulation 
volume, no statistically significant pair bias is found. 
Such poor statistics arise because only a small fraction 
(15%) of total mergers are represented by couples at 
such close separations. Trying to increase statistics by 
considering couples at larger separations is undermined 
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by the fact that the fraction of couples that merge is 
a steeply decreasing function of separation. As already 
mentioned, one cannot consider longer merger timcscalcs, 
as this permits a larger fraction of the halo population 
to have undergone a major merger, thereby washing out 
the signal in the comparison population. We find sim- 
ilar results when looking at the pair bias for mass cuts 
from M > 10 13 /i^Mq to M > 2 x 1O 14 /i~ 1 M , and 
thus we conclude that pair bias cannot be used to reli- 
ably predict merger bias observationally, in simulation, 
or through current analytic treatment. 

7. CONCLUSIONS 

Cluster merger statistics may provide insight into the 
nature of hierarchical structure formation and the mech- 
anisms by which the largest coherent objects in the uni- 
verse form. We use large- volume, high-resolution N-body 
simulations to investigate the utility of close spatial pairs 
of galaxy clusters as proxies for cluster mergers. We char- 
acterize merger statistics through the merger kernel, and 
examine the density dependence of merger efficiency. We 
highlight our conclusions as follows: 

• Close pairs of galaxy clusters at very small sepa- 
rations (< 1 — 2 /i -1 Mpc) can be used to reliably 
predict mergers. However, since these constitute a 
small fraction of the total merger population, close 
pairs are not a reliable proxy for cluster merger 
rates. We quantify this by measuring the efficiency 
and completeness of merger candidates identified 
via close pairs as a function of separation and find 
that their intersection typically occurs at a low 
merger fraction (0.5-0.6). 

• We find that close pairs are even poorer proxies for 
mergers between massive galaxy-sized halos. This 
indicates that galaxy pairs will not provide a reli- 
able proxy for galaxy merger rates at high redshift, 
where most galaxies reside in distinct halos. 

• We note that the failure of close pairs as proxies 
for mergers indicates that determination of merger 
rates from spatial statistics, such as the correla- 
tion function, cannot be trusted outside the physi- 
cal size of a single halo. 

• We examine the mass, redshift, and timescale de- 
pendence of pair mergers, finding that the pair- 
merger hypothesis (that close pairs are proxies for 
mergers) at a given comoving separation is most 
accurate at high redshift, high mass, and over long 
merger timescales. In our best scenarios, the in- 
tersection of merger efficiency and completeness is 
at 75%, i.e. 75% of pairs within a given separa- 
tion merge, and these constitute 75% of all merg- 
ers. We also find a nearly universal relation for 
pair merger efficiency and completeness for differ- 
ent mass halos. This relation begins to break down 
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as we approach massive galaxy-size halos, and is 
compromised by redshift space distortions. 

• Redshift space distortions have a devastating im- 
pact on detecting close galaxy cluster pairs in sur- 
veys; nearly all of the merger candidates identified 
in redshift space do not merge. Although an ex- 
trapolation, we expect these results to be robust 
for galaxy-size halos at high redshift. 

• We present the first fit from simulation to the 
merger kernel — a means to describe halo merger 
rates via the halo mass function (coagulation). 

• The merger kernel exhibits dependence on local 
(~ 17/i _1 Mpc) density. Specifically, halo merging 
in our high-mass regime is more efficient in under- 
dense regions. 

• Pairs at large separations (> 3h~ 1 Mpc) are more 
likely to merge in overdense regions. For pairs at 
small separations, low-mass halos are more likely 
to merge in underdense regions, while high-mass 
halos exhibit no environmental dependence. 

• We sought to use cluster pairs to measure the 
merger bias by using the pair bias as a proxy for 
merger bias. We extended the treatment of previ- 
ous analytic work to include the fact that mergers 
result in mass gain; when computing the bias of 
recently merged halos, the comparison population 
should be a set of halos of the same mass as the 
children, instead of the parents. Close spatial pairs 
that reliably merge are too rare to produce a sta- 
tistically significant measure of merger bias, even 
in a (2.2 ft, _1 Gpc) 3 simulation volume. 

In conclusion, we have shown that close spatial pairs 
of galaxy clusters are of limited value as a probe of over- 
all cluster merger rates. We have determined the merger 
kernel for halo coagulation for the first time from simula- 
tions, finding that a statistical description of halo merg- 
ers is of more promise. Further work is needed to extend 
our parametrization of the merger kernel to lower masses 
and higher redshifts, and to explore whether the environ- 
mental dependence of the merger rate persists in these 
regimes. 

We thank G. Jungman and R. Yan for enlighten- 
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Kamionkowski, and M. White for valuable comments on 
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